Interacting quantum rotors in oxygen-doped germanium 
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We investigate the interaction effect between oxygen impurities in crystalline germanium on the 
basis of a quantum rotor model. The dipolar interaction of nearby oxygen impurities engenders 
non-trivial low-lying excitations, giving rise to anomalous behaviors for oxygen-doped germanium 
(Ge:0) below a few degrees Kelvin. In particular, it is theoretically predicted that Ge:0 samples 
with oxygen-concentration of 10 17_18 cm -3 show (i) power-law specific heats below 0.1 K, and (ii) a 
peculiar hump in dielectric susceptibilities around 1 K. We present an interpretation for the power- 
law specific heats, which is based on the picture of local double-well potentials randomly distributed 
in Ge:0 samples. 

PACS numbers: 65.40.Ba, 65.60.+a, 61.72.Bb, 77.22.Ch 
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I. INTRODUCTION 

Rotational states of impurity atoms in crystalline 
solids affect the essential properties of solids at low tem- 
peratures. This is because the rotational degree of free- 
dom of an impurity yields low-energy excitations that 
manifest themselves such low temperatures. A typical 
example is an interstitial oxygen impurity in covalent 
semiconductors such as Si 1-4 , Ge 5 ~ 7 , and GaAs 8-10 . 
The microscopic structure around oxygen impurities in 
those systems has been thoroughly studied in past. As- 
grown Ge, for instance, is usually contaminated by more 
than 10 oxygen atoms per cubic centimeter with a spa- 
tially random distribution. Oxygen breaks a covalent 
bond between two Ge atoms, establishing a puckered 
Ge-O-Ge segment. At low temperatures, those intersti- 
tial oxygen atoms are quantum-mechanically delocalized 
in the annulus around the original bond center of Ge- 
Ge 11-13 . The rotational degree of freedom of oxygen 
yields low-lying excitations in the far-infrared spectral 
region, which has been observed in phonon spectroscopy 
measurements 11 . The energy spectra deduced from these 
measurements can be accurately reproduced by a quan- 
tum rotor model 13-15 ; Each rotating Ge20 unit within 
this model is mapped onto a quantum rigid rotor rotat- 
ing in a two-dimensional plane that is perpendicular to 
the Ge-Ge axis. 

Since the electronegativity of oxygen is larger than that 
of germanium, puckered segments of Ge-O-Ge carry elec- 
tric dipole moments. In the dilute limit, the interac- 
tion between the dipole moments is negligible because 
the dipolar interaction decreases with the third power of 
distance. However, for actual concentrations of oxygen 
p ~ 10 17 — 10 18 cm -3 , a random distribution of oxygen 
impurities causes the occurrence of clusters of nearby 
dipole moments with a large coupling energy. If the 
coupling is strong, such that it is sufficient to interfere 
with the nearly free rotation of individual Ge2 units, 
the quantum nature of the Ge2 O units that are involved 
in the clusters is completely different from those of the 
isolated Ge20 units. Recently, anomalous dielectric re- 
sponses in coupled dipolar rotors have been reported 16,17 , 



wherein the anomalies originate from peculiar low-lying 
excitations of the interacting dipolar rotors. Therefore, 
it is expected that ensembles of clustered Ge20 units 
play an important role in the low-temperature proper- 
ties of Ge:0; specific heats and dielectric susceptibilities 
are cases in point. 

The present article theoretically investigates the ef- 
fect of the dipolar interaction of rotating Ge20 units on 
the low-temperature properties of Ge:0. It is shown that 
the strong dipolar interaction between clustered Ge2 
units results in peculiar low-energy states, which in turn 
give rise to a power-law temperature dependence in spe- 
cific heats C(T) oc T 5 below 0.1 K under our numeri- 
cal conditions. Further, dipolar transitions between the 
low- lying levels of the interacting Ge20 units engender a 
non-trivial maxima in the dielectric susceptibilities x(T) 
around 1 K. The power-law behavior of specific heats 
is described by a scenario based on the presence of lo- 
cal double-well potentials in crystalline Ge. It should be 
noted that this picture of crystalline solids is exceedingly 
analogous to the theoretical model, referred to as the 
two-level tunneling model, which describes the power-law 
specific heats in amorphous solids 18 . 

The outline of this article is as follows. In Sec. II, 
we introduce a quantum rotor model that accounts for 
the quantum rotation of the Ge20 unit and formulate 
the Hamiltonian for interacting quantum rotors that are 
coupled via dipolar interaction. Section III focuses on the 
low-energy excitations of paired and clustered rotors; the 
distribution function P{e) of the lowest excitation energy 
e is evaluated in line with the theory of the two-level tun- 
neling systems. Section IV shows the calculated results of 
the specific heat C(T) for Ge:0. The power-law behavior 
of C(T) is understood by taking into account the distri- 
bution function P(e) of strongly coupled rotors. Section 
V describes the results of the dielectric susceptibilities 
x(T) of Ge:0. An anomalous hump in xCO around 1 K 
is understood by considering the selection rule for dipo- 
lar transition. Section VI comprises remarks on the rele- 
vance of our findings to other systems. The final section 
concludes the present article. 
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II. SYSTEM 

A. Isolated quantum-rotor model 

We first describe the microscopic structure of an iso- 
lated Ge2 unit, taking no account of interaction be- 
tween them. Figure 1 (a) illustrates the atomic dispo- 
sition of a puckered Ge-O-Ge segment together with six 
Ge atoms next nearest to the interstitial oxygen. At suf- 
ficiently low temperatures, the oxygen atom is quantum- 
mechanically delocalized around the Ge-Ge axis. The 
delocalized oxygen is subjected to repulsive force from 
the six next-nearest neighboring Ge atoms, thus experi- 
encing an azimuthal potential having six-fold rotational 
symmetry 11 . Hereafter, we call it a "hindering poten- 
tial" . 

Recent isotope-shift measurements 13 have revealed 
that the two Ge atoms neighboring O are not fixed, but 
loosely bound around an external axis. As a consequence, 
all three atoms (Ge and O) rotate in phase around the 
external axis as shown in Fig. f (b). The difference of 
the electronegativity between Ge and O causes the oc- 
currence of an electric dipole moment in the direction 
perpendicular to the Ge-Ge axis. This allows us to map 
a rotating Ge20 unit onto a quantum rotor having the 
moment of inertia 11,13 I — f0.6 amu-A 2 and an electric 
dipole moment 19 p = I Debye . 

The Hamiltonian H for an isolated rotor consists of 
two parts: 

H = K + V. (f) 
The first term K accounts for the kinetic energy of the 



rotor associated with its rotation, given by 

h 2 d 2 



K = — 



21 dd 2 ' 



(2) 



The quantity E K = h 2 /(2I) det< :rmines the energy scale 
of the quantum rotation. The second V corresponds to 
a hindering potential, which reads 



V = V cos(66> 



(3) 



with a constant <j>. In the following, we fix Ek — 2.3 K 
and Vo = 1.5Ek in accordance with the values deduced 
from the spectroscopy measurement 11 . The Schrodingcr 
equation for the Hamiltonian (f ) can be solved analyti- 
cally by mapping it onto the Mathieu equation; the de- 
tails of these calculations are given in Refs. 20 and 21 . 



B. Dipolar interaction of quantum rotors 

When quantum rotors get sufficiently close to each 
other, they can no longer be regarded as being isolated, 
thus, interaction between rotors should be taken into ac- 
count. Let us consider two rotors with dipolar moments 
Pi and pj separated by a distance vector R. Expanding 
the interaction potential in terms of 1/R, the lowest-order 
approximation yields the dipolar term 



i (Pi Pj (Pi ■ R)(pj ■ R) 
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i? 3 



R 5 



(4) 



Here e r is the dielectric constant of Ge crystals. Dis- 
carding a numerical factor, the energy scale of dipolar- 
interaction potential (4) is characterized by the quantity 



p- 





(9 



FIG. 1: (a) The atomic coordinate of a puckered Ge20 unit 
together with six Ge atoms next nearest to the interstitial 
oxygen, (b) Microscopic details of a rotating Ge2 unit. The 
whole Ge2 unit rotates in phase around an external axis. 
The difference of the electronegativity between Ge (<5 + ) and 
O (2<5~) causes the occurrence of an electric dipole moment 
p in the direction perpendicular to the rotational axis. 
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For the case Ek ^ J, the rotors are nearly isolated so 
that the interaction potential (4) can be treated as a 
perturbation. When Ek < J, on the other hand, the 
quantum nature of coupled rotors is completely different 
from that of isolated ones because of their strong dipolar 
interaction. 

At finite temperatures, the thermal fluctuation sup- 
presses the dipolar interaction between the rotors. There- 
fore, we introduce the upper cut-off length i£ max defined 
by J(-R m ax) ~ fcfiT, or equivalently, 



Rma,x — 



P 



A-Ke r k B T 



(6) 



Within this cut-off length, the two rotors correlate via 
dipolar interaction. The value of i? max depends on the 
temperature T. For example, T = lOmK results in the 
cut-off length R 

max rsJ 

30A. For actual oxygen concen- 
tration p = I0 17 — I0 18 cm~ 3 , this length is much less 
than the mean separation of adjacent rotors in Ge:0. 




FIG. 2: Definitions of the coordinates of two quantum rotors 
carrying dipole moments pi and p2- The two rotors are sep- 
arated by a vector R and interact with each other via the 
dipolar interaction potential Wi2- 

The length i? max further shortens with increasing tem- 
peratures. Therefore, most rotors may be considered as 
non-interacting rotors above the mK range. However, 
it should be noted that the inhomogcncous distribution 
of rotors makes the formation of clusters consisting of 
two or more rotors. These clustered rotors yield pecu- 
liar low-energy excitations originating from the dipolar 
interaction between rotors, which significantly affect the 
properties of Ge:0 at 1 K and below, as will be shown 
in the subsequent sections. Hereafter, we fix the cut-off 
length -R m ax = 30 A for the temperature range under 
consideration; for simplicity, we do not consider the dis- 
creteness of the lattice of Ge crystals. 

Before closing this section, the spatial correlation of 
impurity distribution should be mentioned. We have as- 
sumed that no correlation exists in the distribution of 
oxygen impurities. Yet this assumption is not always 
correct. When the sample is prepared by pulling from 
the melt, the correlation may be established at tempera- 
tures where oxygen impurities can diffuse. Subsequently 
the correlated distribution of oxygen would be fixed on 
quenching. It is difficult to definitely account for this 
correlation since it depends essentially on the history of 
the sample preparation. Therefore, as a rule, we take a 
random distribution for oxygen with no spatial correla- 
tion, bearing in mind that it is valid for the case of the 
absence of correlations. 



III. LOW-ENERGY EXCITATIONS OF PAIRED 
ROTORS 

In this section, we focus mainly on the quantum char- 
acter of paired rotors (n = 2). This is because the be- 
havior of clustered rotors with n > 3 is essentially under- 
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FIG. 3: Spatial profile of the potential field F in the (61,62) 
plane. The dark regions correspond to negative values of the 
potential F, while the white regions correspond to the pos- 
itive. Two global minima (maxima) are indicated by white 
(black) dashed circles. 

stood by referring to those of paired rotors. Emphasis 
will be placed on the fact that strongly interacting rotors 
give rise to peculiar low-lying excitations; the eigenen- 
ergies of these excitations are well below the energies of 
first-excited states for an isolated rotor. We will see that 
the low-lying excitations of paired rotors markedly con- 
tribute to low-temperature properties of Ge:0 below 1 
K. 



A. Paired-rotor system 

Suppose that two quantum rotors with dipole moment 
Pi(i = 1,2) are positionally separated by the vector R as 
depicted in Fig. 2. The Hamiltonian for paired rotors is 
given by 

H = K 1 +K 2 + V 1 + V 2 + W 12 , (7) 

where the definitions of Ki, Vi, and Wij are the same as 
those in Eqs. (2), (3), and (4), respectively. The explicit 
form of W12 as a function of (6-1,62) is given in Ref. 16 . 

As seen from Eq. (7), the quantum state of paired ro- 
tors is identified with that of single quantum particle 
moving in a potential field 

F(6x ,62) = Vi (6!) + V 2 (62) + W 12 (6i ,6 2 ). (8) 

For later use, we plot the spatial profile of the potential 
F in Fig. 3. The angles (a\,a 2 ,(3) defined in Fig. 2 are 
taken to be (7r/4, — 7r/6, 7r/3), and the interaction energy 
J is taken to be equal to the amplitude of the hindering 
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potential. The dark regions correspond to negative values 
of the potential F, while the white regions correspond to 
the positive. A spatial fluctuation with a period 7r/3 
stems from the hindering term Vi oc cos 60,. We notice 
that two global minima (maxima) appear at the anti- 
parallel (parallel) dipolar configurations as indicated by 
a white (black) dashed line. These global minima and 
maxima indicate that two dipoles prefer to adopting an 
anti-parallel configuration. The positions of two minima 
(and maxima) generally deviate from the lines \9 2 — 9\ = 
and \9 2 — #i I = k due to the effect of the finite angles 
(oei,a 2 , f3). For any set of angles {ot\,a 2 , f3), however, the 
relative difference between the position of the minimum 
and the maximum are invariant with A0i = A02 = n. 

B. Eigenenergies for paired rotors 

The Schrodinger equation H^(9 1 ,9 2 ) = EV(6 1 ,6 2 ) for 
the Hamiltonian (7) is solved by taking the solution in 
the form 

*(0i,02) = ^Q 1 ,; 2 exp{-i(/ 1 1 + Z 2 2 )}, (9) 

where (i = 1,2) takes an integer value h = 
0, ±1, ±2, ±3, • • • . The coefficients Ci lt i 2 are obtained nu- 
merically by diagonalizing the finite matrix of H in the 
subset of basis states with |Zi < l\ and \l 2 \ < l 2 . The 
cut-off values Zf and l 2 are increased until the considered 
eigenvalues converge within the desired accuracy. For ac- 
tual calculations, we set Zf = l 2 = 30 in order to obtain 
precise low-energy states. 

Figure 4 shows the eigenenergy spectra for paired ro- 
tors as a function of the dipolar interaction energy J. 
Angles (a\, a 2 , f3) were chosen to be (ir/4, —ir/6, 7r/3) as 
an example. When J = 0, the eigenenergies are almost 
degenerate at E/E K = l\ + l\ = 0, 1, 2, 4, 5, 8, • • • . An 
increase in J leads to a split off of the degeneracies at 
J = 0, followed by a spread along the ordinate. Note 
that the energy difference between the ground state and 
the first-excited state decreases monotonously with in- 
creasing J. For large J > Ek, these two eigenlevels form 
a doublet with very small energy splitting. Besides, a 
part of higher eigenlevels also form doublets at J > Ek 
as shown in Fig. 4. As a consequence, a few doublets lo- 
cate within the low energy region [Eg, Eq + Ek], where 
Eg is the eigenenergy of the ground state. Those dou- 
blets locating at the above energy region are what we call 
"very-low-energy excitations" of paired rotors. We notice 
that the very-low-energy excitations of paired rotors are 
generated only when the rotors are strongly coupled. 

The occurrence of the doublets for J > Ek is at- 
tributed to the spatial profile of the potential field 
F(9i,6 2 ) defined in Eq. (8). For strongly paired rotors 
satisfying the condition J > Ek, the depth of the two 
minima in the 0i-0 2 plane become so large that the field 
F eventually forms a double- well potential (See Fig. 3). 
As a result, the behavior of paired rotors is described 
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FIG. 4: The eigenenergy spectrum for paired rotors as a func- 
tion of the interaction energy J. The relative angles of rotors 
are taken to be (ai,Q2,/3) = (w/4, — 7r/6, 7r/3), as an example. 

by a quantum particle moving in a double-well potential 
whose barrier height is the order of the interaction en- 
ergy J. When the kinetic energy of the particle is much 
smaller than the barrier height, eigenfunctions at low en- 
ergies are described by a superposition of two wavefunc- 
tions localized in the respective wells. A slight overlap 
of the localized wavefunctions results in a doublet with 
small energy splitting. By raising the interaction energy 
J, therefore, the barrier height becomes so large that 
the energy splitting monotonously decreases as shown in 
Fig. 4. 



C. Low-energy excitations in clustered rotors 

Doublets locating at very-low energies also occur in 
clustered rotors with larger sizes n > 3. Figure 5 gives 
the energy spectra of the clustered rotors; twenty clus- 
ters of various sizes and various spatial configurations 
were randomly chosen. For any clusters, the ground 
state energies are fixed at zero for the sake of illustration. 
Most eigenlevels are located densely at specific energies 
E/Ek — 0, 1, 2, • • • , i.e., at eigenenergies of an ensemble 
of isolated rotors. This indicates that the rotor's interac- 
tion within individual clusters is not significantly large. 
Nevertheless, a limited number of eigenlevels form dou- 
blets showing energy splitting much smaller than Ek- 
These doublets originate from the small separation R of 
the rotors involved in each cluster. Namely, very-low- 
energy excitations in clustered rotors can be found in 
clusters of any size because of a statistical fluctuation of 
the rotor's spatial distribution. 
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FIG. 5: Eigenenergy distributions for clustered rotors. 
Twenty clusters were randomly chosen with various sizes and 
various spatial configurations. 



Recall that the value of the kinetic energy Ek for an 
actual Ge20 unit equals 2.3 K. This implies that only 
the very-low-energy excitations satisfying the condition 
E/Ek < 0.5 are relevant to the properties of Ge:0 below 
1 K. Further, when considering the properties of Ge:0 
far below 1 K, the presence of doublets with the lowest 
energies would be especially important. In the next sub- 
section, we discuss the energy distribution of such the 
lowest-energy excitations, which would contribute to the 
specific heat of Ge:0 far below 1 K. 



D. Two-level tunneling systems in Ge:0 

As mentioned before, the occurrence of peculiar very- 
low-energy excitations stems from the potential field 
F(9i,02) defined in Eq. (8). For strongly coupled ro- 
tors, the field F forms a double- well potential in the d\- 
#2 plane, and the quantum behavior of the rotors is de- 
scribed by a quantum particle moving in the double-well 
potential. Thereby, strongly coupled rotors yields locally 
a quantum tunneling system. Those tunneling systems 
are randomly distributed in Ge:0 due to an inhomoge- 
neous distribution of oxygen impurities. The height of 
the potential barrier U separating the two potential min- 
ima depends on the distance R between rotors as well 
as the internal configuration of interacting rotors. Dis- 
carding numerical factors, the energy scale of the barrier 
height U is estimated of the order of the interaction en- 
ergy J. 

The existence of tunneling systems in crystalline solids 
has also been demonstrated in alkali halides with sub- 
stitutional impurities such as KCl:Li, NaCLOH, and 
KBr:CN. 22 In these solids, tunneling occurs due to a sub- 
stitutional single defect sitting in sites offered by the host 
lattice. Each defect has an electric or elastic dipole mo- 
ment. The possible orientation of the dipole is deter- 



mined by the local potential minima in crystalline en- 
vironment surrounding the defect. As a result, the de- 
fect atom can tunnel through several equivalent potential 
minima. 

Regarding KBr:CN, Solf and Klein 23 have theoretically 
investigated the energy spectra of the system with paired 
tunneling dipoles. Polar molecules of CN~ dissolved in 
KBr form interacting dipoles oriented in eight directions. 
For 340 ppm of CN~ in KBr, they showed that the dis- 
tribution function P(s) of the lowest excitation energy 
e for paired elastic dipoles changes from P(e) ~ const, 
for e/k B < 100 mK to P(e) cx e" 1 / 2 for e/k B < 10 
mK. 23 This yields a power-law temperature dependence 
of the specific heats as C(T) oc T at 10 mK< T <1 K 
and C(T) oc T 1 ' 2 at 1 mK < T < 10 mK. The point 
is that paired tunneling impurities engender novel low- 
energy excitations with small energy splitting e. This is 
actually realized for interacting quantum rotors in Ge:0 
as mentioned earlier, though the explicit form of P(e) 
remain unclear. The microscopic character of interacting 
constituents in Ge:0 (weakly hindered quantum rotors) 
is evidently distinct from that in alkali halides (n-oriented 
tunneling dipole). Nevertheless, the difference is not in- 
fluential to the emergence of the low-energy excitations 
with small energy splittings. 

Hereafter, we discuss the energy dependence of the dis- 
tribution function P(e) of the lowest energy excitations 
e for strongly paired rotors in line with the theory of 
two-level tunneling systems. 24 ' 25 The height of the tun- 
neling barrier U between two potential minima scales as 
U ~ J cx i? -3 . Thereby the WKB approximation 26 al- 
lows us to evaluate the tunneling amplitude A p express- 
ing the coupling between two localized wavefunctions as 



A p - 7 exp (-aU 1/2 ^j - 7 exp (-aiT 3/2 



(10) 



The distribution function P(e) is given by the following 
three dimensional integral: 

P(e) ~ J dR5{e- A P (R)} . (11) 



Substituting Eq. (10) into Eq. (11), we have 
P{e) cx 



•{log(7/<r)} 



(12) 



For strongly coupled rotors, the energy splitting e is much 
smaller than the parameter 7. When e <C 7, the function 
P{e) given in Eq. (12) rapidly increases with decreasing 
e. Its asymptotic behavior is approximated by the power 
law of P(e) cx £~ q , where the exponent a is of the order 
of unity or less. 

The ensemble of the two-level tunneling systems that 
obey the distribution function P(e) contributes to the 
specific heat C{T) by the following formula: 24,25 



(13) 



C(T)= / deP(e)Co(e,T), 
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and 



C (e,T) oc 



k B T 2 



sech 



2k R T 



Hence, the asymptotic behavior of P(e) oc e~ a results 
in the power-law temperature dependence of the specific 
heats C(T) oc T 1_a . This implies that the specific heats 
for Ge:0, as in the case of KBr:CN, should exhibit power- 
law behavior at sufficiently low temperatures. 

It should be noted that the exponent of the predicted 
specific heat C(T) oc T 1- " for Ge:0 does not need to co- 
incide with the exponent of the power-law specific heat 
for KBr:CN. This is due to the difference of the micro- 
scopic structure of tunneling constituents between in the 
two systems. For KBr:CN, it has been assumed 23 that 
the local potential minima that determines the possible 
orientation of dipoles are invariant to the increase in the 
interaction energy J. This assumption has led to the 
consequence that the lowest excitation energies e are pro- 
portional to J -1 and J~ 2 , within the second- and third- 
order perturbation theory, respectively. As a result, the 
e-dependence of the density of states P{e) can be solved 
analytically, which leads to a linear or T 1 / 2 specific heats 
depending on the temperature range. In contrast, for 
interacting quantum rotors in Ge:0, the energitically- 
preferable orientation of rotors complicatedly depend on 
both J and the spatial configuration of interacting ro- 
tors, namely, by the angles (ai, a%, [3). Hence, the lowest 
excitation energies e can not be expressed as a simple 
function of J and (a\, a.2, /?), which prevent us from ob- 
taining analytically an explicit form of P(e) and C(T). In 
order to clarify the T-dependence of the low-temperature 
specific heats C(T) in Ge:0, therefore, it is crucial to nu- 
merically simulate C{T) directly. 



Here V represents the volume of the system to be consid- 
ered. Equation (16) amounts to considering all clusters 
(14) with different sizes and with different spatial configura- 
tions. This thereby gives the specific heat for a Ge:0 
sample, namely, for a particular realization of the rotor's 
spatial distribution. Since this distribution strongly de- 
pends on the samples used, we have to take sample av- 
erages in order to obtain the specific heat for a typical 
Ge:0 sample. 

To perform the above procedure, we have numerically 
simulated more than 100 samples of Ge:0; in each sam- 
ple, 10 6 rotors are inhomogeneously distributed in a cubic 
system, then are factorized into many n-sized clusters. 
The volume of the system is defined by 10 6 //?, and the 
lower cut-off i? m i n = 3 A is introduced in order to pre- 
vent two (or more) oxygen impurities from getting closer 
than the bond length of Ge-O-Ge. 



IV. SPECIFIC HEATS 
A. Method 

Now let us consider the specific heat for a Ge:0 sam- 
ple having an impurity concentration p. This sample 
contains a large number of quantum rotors; most rotors 
are isolated because of large separations, while a part of 
neighboring rotors establish clusters having various sizes 
n. The specific heat C c \ s for an individual cluster is cal- 
culated straightforwardly through its partition function 
Zeis = tr(e- fJH ) and internal energy U c \ s = -^logZ cls , 



a 



els 



dU c \ s 
dT ■ 



When calculating specific heats for a whole Ge:0 sample, 
on the other hand, we have to average the contributions 
of all clusters involved in the system: 



els 



B. Numerical results 

Figure 6(a) shows the numerical results of specific heat 
C(T) for various concentrations of rotors p: (i) p = 10 18 
cm~ 3 , (ii) 3 x 10 17 cm~ 3 , and (iii) 10 17 cm" 3 . The dipole 
moment p of quantum rotors is set to be 1 Debye as es- 
timated before. We also plotted the Debye specific heat 
C(T) = aT 3 with a = 5.83 x 10~ 2 g" 1 K" 4 . The 
value of a is evaluated based on the experimental data ob- 
served at the temperature range 2.8-100 K 27 . In Fig. 6 
(a), broad peaks at around T = 1 K commonly appear 
for all rotor concentrations. These peaks are caused by 
the fact that most rotors contained in the samples are 
independent of each other. For isolated rotors, the level- 
splitting between the ground state and the first-excited 
ones is equal to Ek = 2.3 K as seen in Fig. 4. Thus the 
ensemble of isolated rotors contributes to specific heats 
as the Schottky peaks at around Ek/2 ~ 1 K. 

Of particular interest is the temperature dependence 
of specific heats C(T) below the Schottky peaks. At 
temperatures T < 0.1 K, the calculated results obey a 
power-law temperature dependence of the form oc T - 5 . 
The exponent 0.5 seems to be independent of the im- 
purity concentration p, whereas the magnitude of C(T) 
shifts downward with decreasing p. Furthermore, the 
magnitudes of C(T) for temperatures T < 0.1 K are 
proportional to p 2 , while those for T > 0.1 K (around 
the Shottky peak) show a trivial linear dependence on p. 
Interestingly, the same p 2 -dependence of the power-law 
specific heats has been pointed out in KBr:CN 23 . The 
physical origin of these peculiar behaviors in C(T) will 
(15) be addressed soon below. 

We have also clarified the dependence of C(T) on the 
definition of the value of p for individual rotors. Fig- 
ure 6 (b) shows the resulting C(T) for p = 1 Debye 
(solid line), 0.75 Debye (dashed-dotted), and 0.5 Debye 
(dashed). Decreasing the value of p reduces the magni- 
fy) tudc of C(T), and results in a deviation of C(T) from 
the power-law temperature dependence. Nevertheless, 
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FIG. 6: (a) Calculated results of specific heats C(T) for Ge:0. The concentration of rotors p is varied as follows: (i) p = 10 18 
cm ~ 3 , (ii) 3 x 10 17 cm" 3 , (iii) p = 10 17 cm , with the dipole moment fixed at p = 1 Debye. The dashed-dotted line follows 
the Debye specific heat for pure Ge crystals deduced from Ref. 27 . (b) Results of specific heats C(T) for p — I0 18 cm~ 3 with 
varying the definition of p; (i) p = 1 Debye, (ii) 0.75 Debye, and (iii) 0.5 Debye. 



the considerably excess specific heat below 0.1K remains 
observed for all ps, which is a manifestation of the dipole 
interaction of rotors consisting of clusters. It should be 
emphasized that, for all numerical conditions in Figs. 6 
(a) and 6 (b), the results of C{T) below 0.1 K well over- 
come the Debye specific heat. This demonstrates that 
our theoretical results for specific heats, both the Schot- 
tky peaks and the power-law behavior, can be observed 
in experiments. 



C. Two anomalies in the excess specific heats 
below 0.1 K 

In Figs. 6 (a), we have found two striking features in 
the excess specific heats below 0.1 K. The one is the 
power-law temperature dependence of C(T) oc T 5 , and 
the other is the p 2 -dependence of the magnitude of C(T). 
To reveal the physical origin of the two features, we have 
clarified the contribution of the ensembles of clustered 
rotors of respective sizes. 



1. Power-law specific heats 

In Fig. 7, the result of C{T) for p = 10 18 cm" 3 and 
p = 1 Debye is replotted by a thick line. In addition, we 
exhibit four components consisting of C(T) for a com- 
plete Ge:0 system; solid thin lines show the contribution 
of the ensemble of coupled rotors (n = 2), and dashed 



lines show that of clustered rotors (n > 3), as denoted 
in the figure. The two components, indicated by the 
term "weak" , correspond to the ensemble of clusters es- 
tablished by weakly-interacting rotors; within those indi- 
vidual clusters, all rotor separations are larger than the 
characteristic value R c = 10A (See Fig. 8 (a)). Because 
of the large separations, the eigenenergies of those clus- 
ters are nearly equal to the eigenenergies of isolated ro- 
tors, thus giving rise to the Schottky behavior appearing 
in Fig. 7. On the other hand, the other two components, 
indicated by "strong" , represent the contribution of the 
ensemble of clusters including strongly-interacting rotors; 
within those clusters, two or more rotors are closer than 
R c (See Fig. 8 (b)). As clearly seen from Fig. 7, the excess 
specific heats below 0.1 K originate from the ensemble of 
strongly coupled rotors (n = 2) with separations R < R c . 
Consequently, we conclude that the power-law behavior 
of C(T) cx T - 5 for a Ge:0 sample originate from the 
very-low-energy excitations of strongly coupled rotors. 

Our conclusion mentioned above is justified by eval- 
uating the behavior of P(e) that directly relates to the 
low-temperature specific heats through Eq. (14). Fig- 
ure 9 shows the distribution function P{e) of the energy 
splitting e between the ground state and the first-excited 
ones for strongly coupled rotors; the definition of p is set 
to be p = 1 Debye for solid line, 0.75 Debye for dashed 
line, and 0.5 Debye for dotted line. All numerical re- 
sults are normalized as J P(e)de = 1. The definitions 
-Rmin = 3A and R c = 10A give a lower limit for the en- 
ergy splitting £ m ; n and an upper one e max , respectively. 
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Temferature T[^ 

FIG. 7: The result of C(T) for p = 10 18 cm" 3 and p = 1 
Debye (thick line) is replotted, together with four compo- 
nents stemming from ensembles of paired (n = 2) and clus- 
tered (n = 3) rotors. The meanings of the term "weak" and 
"strong" are given in text. 
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FIG. 8: Schematic illustrations of rotor clusters with n = 
3: (a) A cluster formed by weakly interacting rotors. All 
separations Rij exceed the characteristic length R c — 10A. 
(b) A cluster established by strongly interacting rotors. The 
separation R12 is smaller than R c . 



2. p 2 -dependence of excess specific heats 



Numerical results of P(e) show a maximum at around 
1.5-2.0 K followed by a rapid decrease with decreasing 
e. Interestingly, for a lower e than 0.1 K, the function 
P(e) (except for dotted line) obey the power-law form 
P{e) oc s~ a within a range over one order (See the inset 
of Fig. 9). This resulting power-law is consistent with the 
theoretical consequence given in Sec. III-D. In the case 
of p = 1 Debye, in particular, the exponent a takes 0.5; 
this reproduces the power-law specific heat C(T) oc T 5 
for sufficiently low temperatures T <C £ m ax/fc_B, which is 
clearly demonstrated in Fig. 6 (a). 

We have numerically confirmed that the exponent 0.5 
in the power-law form of P(e) oc e -0 - 5 is invariant to the 
change in the rotor concentration p under the condition 
p = 1 Debye. On the other hand, for smaller p, the 
exponent a of P(e) oc e~ a for e/k B <C 0.1 K gradually 
decreases, and finally P(e) becomes a constant as shown 
by the dotted line (for p = 0.5 Debye) in the inset of 
Fig. 9. This results from the fact that, when p decreases, 
the lower bound e m in increases because of the reduction of 
the interaction energy J. Eventually £ m ; n gets larger than 
the energy range where the function P(e) clearly exhibits 
a power-law behavior P{e) oc e~ a with positive a. As a 
result, P(e) seems to be constant for the condition of 
p = 0.5 Debye. In addition, the energy range where P(e) 
is constant is rather narrow; this is the reason why the 
specific heats for the condition of p = 0.5 Debye exhibit a 
power-law temperature dependence only within a rather 
narrow temperature range as shown in Fig. 6 (b). 



Next we turn to the p 2 -dependcnce of the magnitude 
of C(T) below 0.1 K. The p 2 -depcndence of the low- 
temperature specific heats has also been pointed out for 
KBr:CN as a consequence of the virial expansion of the 
free energy up to p 2 23 . In the following, we derive 
the /9 2 -dependence of C(T) straightforwardly by consid- 
ering the number of strongly coupled rotors N p in the 
system. The quantity N p can be expressed by the re- 
lation N p — N x u. Here, N is the total number of 
rotors in the system, which is apparently proportional to 
the rotor concentration p. The quantity u represents the 
possibility that a given rotor finds another rotor with a 
separation R < R c . Thereby u reads 



u= [ P(R)dR, (17) 

where P(R) is the distribution function for the separation 
R. 

The explicit form of the function P(R)dR is obtained 
as follows 28 . P{R) represents the probability of finding 
no rotor within a sphere of radius R and one rotor in the 
spherical shell of the thickness dR, when one rotor is lo- 
cated at the origin. In order to calculate P(P), we divide 
the sphere of radius R into m shells of thickness R/m. 
Since the rotors are randomly distributed, the probability 
of not finding a rotor in a sphere of radius R is expressed 
by the probabilities of none of these m shells contain- 
ing rotors. Thus, we can express the probability density 
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FIG. 9: Distribution functions of the energy splitting e be- 
tween the ground state and the first-excited one. The defini- 
tion of p is set to be p = 1 Debye for solid line, 0.75 Debye 
for dashed line, and 0.5 Debye for dotted line. 



P(R) for large m 
P(R)dR = 



1 4 ( R 

1 - T^P — 

3 \m 



rrm-l 

X LL j=o 



1-Anp[^) ■ * 
ml m 



x 4irpR 2 dR. (18) 



Taking the limit m — > oo, the probability of finding a 
nearest- neighbor rotor between R and R + dR is given by 



P(R)dR = Att P R z cxp \--npR J dR. (19) 

The probability density P(R) has its maximum at Rq = 
(27T/9) -1 / 3 . Since R c = 10 A is much smaller than R , the 
function P{R) can be approximated by P(R) ~ A-xpR 2 . 
Substituting it into Eq. (17), we obtain the relation u oc 
p. As a result, the number of strongly coupled rotors N p 
is proportional to p 2 . This naturally leads to the peculiar 
/? 2 -dependence of the magnitude of C(T) below 0.1 K, as 
demonstrated in Fig. 6 (a). 

Let us summarize the findings for low-temperature spe- 
cific heats in Ge:0 systems. Numerical simulations have 
suggested the three anomalous behaviors in specific heats 
of Ge:0 systems: (i) the Schottky peaks at around 1 K, 
(ii) the power-law temperature dependence below 0.1 K, 
and (iii) the p 2 -dependence of the magnitude of C(T) 
below 0.1 K. The first can be attributed to the ensem- 
ble of isolated rotors in the system, while the latter two 
stem from the presence of strongly coupled rotors with 
small separations R < R c = 10A. Furthermore, the phys- 
ical origin of the latter two anomalies has been quanti- 
tatively understood through the distribution function of 
the energy splitting P(e) as well as that of separation 
P{R)- It is expected that our predictions will shed light 



on experimental research regarding the low-temperature 
properties of oxygen-doped semiconductors. 



V. DIELECTRIC SUSCEPTIBILITY 

We now focus on the dielectric response for Ge:0 at low 
temperatures. The quantum tunneling dipolcs in alkali 
halides crystals are known to significantly contribute to 
the dielectric susceptibility at low temperatures. In a se- 
ries of papers 29 , Klein has theoretically investigated var- 
ious dynamic and static quantities in terms of the virial 
expansion, suggesting that the low-T dielectric suscepti- 
bility contains a — logT term. The experimental 
results of x(^) f° r KCbLi have exhibited an interesting 
dependence on the impurity concentration p and tem- 
perature T. At p of a few ppm, the susceptibility x(^) 
monotonically increases with decreasing T, and then ap- 
proaches a constant value for T — > 0. On the other hand, 
for the samples of 60 ppm or more, xCH is n0 longer 
a monotonic function of T. The susceptibility x(T) dis- 
plays a maximum at a temperature of approximately 300 
mK, and the magnitude of the maximum is pronounced 
by an increase in the impurity concentration p. Further, 
an additional increase in p of 1100 ppm surprisingly re- 
sults in a decreasing susceptibility, which is a manifesta- 
tion of the nonlinearity of the low-T susceptibility xCH 
on p. In Refs. 30 ' 31 , the occurrence of the maximum in 
x(T) as well as the abovementioned non-linearity of x(T) 
on p is described by the theory based on Mori's reduction 
method 32 ' 33 . 

In this section, we examine the effect of the presence 
of the interacting quantum rotors on the dielectric sus- 
ceptibility x(T) in Ge:0 for an oxygen concentration 
p = 10 17 - 10 18 cm" 3 . The maximum of x(T) at 1 K, 
similar to that observed in KCl:Li, is also observed for 
Ge:0. Interestingly, the origin of the maximum in x{T) 
is directly attributed to the presence of the low-lying ex- 
citations of the interacting quantum rotors, which is dif- 
ferent from the interpretation of the maximum of x(T) 
for KCl:Li, given in Refs. 30 - 31 . 



A. Method 

The linear dielectric function for an individual clus- 
ter is described by the time-dependent response function 
with respect to an external field, 

x?i(t f) = PAnm n (20) 

Here, the angular brackets indicate that the thermal av- 
erage (•••) = Z~ x tr (• • • e~^ ff ) is taken. The operator 
P M (t) represents the ^-component of the total dipole mo- 
ment for the cluster, and 9(£) is the Heviside step func- 
tion. Since experiments provide quantities depending on 
frequency rather than on time, we take the Fourier trans- 
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FIG. 10: Schematic illustration of allowed dipolar transitions: 
(a) the case for an isolated rotor, and (b) single paired ro- 
tors. Vertical lines indicate the allowed dipolar transitions; 
a dashed line refers to the allowed transition involving the 
ground state, and a solid line refers to the remainder. 



form of Eq. (20), 



eZ 



3,1*3 

E4 — Ei 



{Ei-Eif-{fujf 



cxp 



Here, \Ej) means the eigenstates of the clustered rotors 
belonging to the eigenenergy Ej . We focus on the dielec- 
tric response of Ge:Oi at frequencies u lower than the 
kHz range; in this case, the energy difference \Ej — Ei\ 
for arbitrary eigenstates is larger than fko by many or- 
ders of magnitudes. Thus in the sequel we take the zero- 
frequency limit of Eq. (21). 

For calculating the susceptibility x^v for a whole Ge:Oi 
sample, we have summarized the contribution of all clus- 
ters, 



(22) 



els 



and performed the sample average as in the cases for 
calculating the specific heat. In addition, the isotropy 
of Ge:Oi allows us to consider the quantity \ = (x xx + 
Xyy+Xzz)/3. Consequently, the objective is the following 
quantity, 



X(T) = £x cls (T) 
X ds (T) 

with arbitral component P M . 



els 

2 ^ |(£ 3 |P M |^)| 2 
3eZ^ Ej — Ei 6XP 

3-1*3 J 



3l 



(23) 



24) 



B. Contribution of isolated and paired rotors 

For later use, we examine the dielectric properties of a 
isolated quantum rotor and that of paired rotors. Figure 
10 displays the low-lying energy levels of (a) an isolated 
rotor, and (b) weakly coupled rotors. Vertical lines in- 
dicate the allowed dipolar transitions; dashed lines refer 
to the allowed transition involving the ground state, and 
solid lines to the remainder. Dipolar transitions involving 
higher levels than E/Ek = 4 are not shown here, since 
those higher levels are only negligibly excited below a few 
degrees Kelvin. 

For an isolated rotor, only two dipolar transitions 
contribute to the dielectric susceptibility x(^) below a 
few Kelvin. On the other hand, for single paired ro- 
tors, many eigenlevels lower than AEk contribute to the 
dielectric susceptibility as shown in Fig. 10 (b); rings 
bundling four eigenlevels indicate that the correspond- 
ing vertical line connects all of the four eigenlevels bun- 
dled. As a consequence, 24 dipolar transitions are ex- 
hibited in Fig. 10 (b). Particularly important is the oc- 
currence of the allowed transition at around E ~ Ek 
connecting nearly-degenerate eigenlevels. The small en- 
ergy differences \Ej — Ei \ of two nearly-degenerate eigen- 
states \Ej) and \E{) magnify the absolute value of the 

I (E ■ I _P \Ei ) 1 2 

fraction e-Ei a PP earm g m Eq. (24), then leading 
anomalous maximum in x(^) a t the temperature 



Ex/ks ■ We emphasize that such almost degen- 



to 

T > 

erate levels at E ~ Ek emerge only for weakly paired 
(and clustered) rotors. Namely, the remarkable contri- 
bution of the dipolar transitions associated with nearly- 
degenerate levels occurs only for interacting rotors. This 
implies that the presence of paired and clustered rotors 
affects the dielectric properties of Ge:0 at low tempera- 
tures. 

We have confirmed that the maximum in x(T) appears 
also for clustered rotors with sizes n > 3. Note that, 
however, the magnitude of x(T) for clustered rotors is 
generally reduced with increasing the cluster size. This 
is because the many-body interaction of rotors with ran- 
dom configurations prevents them from orienting in the 
direction of an electric field. Indeed, the magnitude of 
x(T) for clusters n — 5 become less than that for paired 
rotors by two orders or more. Furthermore, the possibil- 
ity of the occurrence of clusters with large sizes (n > 6) is 
extremely small. These facts allow us to ignore the pres- 
ence of clusters with n > 6 for considering the dielectric 
susceptibility x{T) in Ge:0 systems. 



C. Dielectric responses for Ge:Oi 

Figure 11 (a) shows the calculated results of the dc di- 
electric susceptibility x(^) f° r Ge:0 samples; (i) p = 10 18 
cm~ 3 , (ii) 3 x 10 17 cm~ 3 , and (iii) 10 17 cm~ 3 with fixing 
the dipole moment p — 1 Debye. The solid lines show 
the total susceptibility x(T), while the dashed ones show 
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FIG. 11: (a) The dc dielectric susceptibility x(^) of the same 
four conditions shown in Fig. 6. Solid lines show the total 
susceptibility x(T), and dashed lines show the contribution 
of the ensemble of isolated rotors. Shallow maxima at around 
T — 1.0 K appear in all conditions, (b) The differences A%(T) 
between the solid lines and the dashed lines exhibited in the 
figure (a). 



the contribution of the ensemble of isolated rotors; the 
latter can be analytically obtained from Eq. (24), since 
the eigenenergics and their eigenfunctions for isolated ro- 
tors are analytically calculated. The magnitude of x(T) 
depends linearly on the concentration of oxygen p and 
squarely on the dipole moment p, which is simply under- 
stood from Eqs. (23) and (24). We find that an increase 
in p gradually enhances the maximum of x{T) at about 1 
K, which comes from an increase in the number of paired 
and clustered rotors in the Ge:0 sample 34 . 

To aid in illustration, Fig. 11 (b) exhibits the differ- 
ences A%(T) between the solid lines and the dashed lines 
displayed in Fig. 11 (a). The maximum in A%(T) is lo- 
cated at about 1 K independent of the oxygen concen- 
tration. Notably, when we compare the data for cases 
(i) and (ii), we find that a change in oxygen concentra- 
tion p by a factor of three enhances the quantity Ax 
by almost an order. Namely, the magnitude of A% de- 
pends almost squarely on the oxygen concentration p. 
This p 2 -dependence of Ax can be explained using the 
same scenario as that for the power-law specific heats 
(See Sec. IV-C). Consequently, the p 2 -dependence of the 
quantity A%(T) as well as the maximum in A\(T) at 
about 1 K directly reflect the presence of paired and clus- 
tered rotors in Ge:0 systems. We emphasize that the 
non-trivial /^-dependence of Ax as well as its maxima 
can be measured experimentally, since the contribution 
of isolated rotors, depicted by the dashed line in Fig. 11 
(a), has an analytic expression deduced from Eq. (24). 

It is worth noting that our results of x(T), shown in 
Fig. 11 (a), are quite similar to the experimental data 



of x(T) for KCl:Li observed at temperatures between 10 
mK and 10 K 30 ' 31 . The experiments were performed 
with Li concentrations ranging from 4 ppm to 1100 ppm, 
revealing the maximum in x(T) between 100 mK and 
500 mK. Further, the peak becomes more pronounced 
with increasing Li concentrations, which is quite similar 
to the case of Ge:0 as shown in Fig. 11. However, de- 
spite the similarity, the physical origin of the maximum 
in x(T) predicted for Ge:0 is different from that observed 
for KCl:Li 30 < 31 . In the latter case, the peak in x(T) is 
identified with the relaxation peak that indicates the col- 
lective motion of the tunneling dipoles 32,33 . On the other 
hand, the maximum in x{T) for Ge:0 systems is not re- 
lated to the relaxation process. The maximum originates 
from the dipolar transition of the nearly degenerate lev- 
els for weakly coupled rotors. In other words, for the 
interacting dipoles, the maximum in x(^) possibly oc- 
curs without the relaxation process. The presence of the 
allowed dipolar transition of the nearly degenerate levels 
is of importance. Consequently the maximum in x{T) 
for Ge:0 that we have observed is different from the re- 
laxation peak observed in KCl:Li. When the oxygen con- 
centration is further increased, the collective motion of 
the rotors would become relevant. In this case, the re- 
laxation effect might manifest itself in the temperature 
dependence of x(T), which is similar to that in KCl:Li. 
Further discussion are required to clarify this point. 



VI. DISCUSSIONS 

In this section, we focus on the relevance of our find- 
ings to other physical systems that also exhibit low- 
temperature anomalies similar to those in Ge:0. 



A. The dipole gap 

It should be remarked the relevance of our model to 
the dipolar gap theory 36 ' 37 . It is known that the long- 
range interaction of constituents in solids generally leads 
to a remarkably reduction of the density of low-energy 
excitations. This reduction stems from the requirements 
of the ground state stability relative to two-particle ex- 
citations and more complex ones for the system of local 
interacting centers. In the case of dipolar interaction, the 
density of states P{e) for low energies tends to be zero 
logarithmically as 36 P{e) cx 1/ log(7/e), where 7 is an ap- 
propriate constant. This soft gap, called the dipole gap, 
manifests itself in a singularity of the specific heats as 
C(T) cx T/log(7/fc B T). Kirkpatrick and Varma 38 have 
suggested the result C(T) oc T 3//2 in terms of the Monte 
Carlo simulations. These conclusions disagree with our 
result of C(T) and P{e) for coupled dipolar rotors. 

We believe that the above disagreement in the behav- 
ior of C(T) is due to the following reasons. The first has 
to do with rather dilute concentrations of rotors in Ge:0 
systems. The occurrence of the dipolar gap requires a 
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large number of dipole moments interacting simultane- 
ously. This is hardly realized in Ge:0 systems, since in 
crystalline Ge there exists the upper limit of solubility 
of oxygen 7 . Namely, a high concentration of rotors suffi- 
cient to realize the dipole gap cannot be achieved as far 
as Ge:0 systems are concerned. 

The second reason is the difference of the nature be- 
tween classical dipole moments and quantum dipolar ro- 
tors. The derivation of the dipole gap in Ref. 36 was based 
on the classical motion of charged particles. Differing 
from that, we have considered the quantum character 
of rotating Ge20 units carrying dipole moments; in ad- 
dition, the rotation of the GC2O units is restricted to 
within a two-dimensional plane perpendicular to the Ge- 
Ge axis. Thus the argument based on the coupled rotor 
model surely yields the density of states P(E) and the 
T-dependence of C(T) distinct from those in Ref. 36 . In- 
dependent of the relevance for Ge:0, nevertheless, it is 
interesting whether or not the dipolar gap occurs in the 
quantum dipolar rotors in sufficiently dense concentra- 
tions. The scaling scheme for disordered spin systems 39 
may help to address this issue. 

B. Two- level tunneling states in amorphous solids 

It is known that the power-law specific heats at low 
temperatures have been observed in various kinds of 
amorphous materials 40 . In those systems, the power-law 
specific heats have been explained by a phenomenological 
model based on two-level tunneling states (TLSs) with an 
assumed constant density of states 24,25 . The presence of 
TLSs in various amorphous solids has been confirmed by 
many experiments 18 . There is, however, as yet no micro- 
scopic derivation for the tunneling states or for the con- 
stant density of states. In addition, the theoretical result 
C(T) oc T differs somewhat from the observed behav- 
ior C(T) cx T 1+s with small deviations S depending on 
materials. Recently found evidences for dipolc-dipole in- 
teractions between TLSs draw attention to the necessity 
of modifying or generalizing the simple TLS picture 41-44 . 

Of particular interest is that the power-law specific 
heat in Ge:0 can be interpreted in line with the sce- 
nario analogous to that of the tunneling model. Empha- 
sis should be made on, nevertheless, the two crucial dif- 
ferences between the coupled-rotor model and the stan- 
dard TLS model as follows: (i) The coupled-rotor model 
gives the clear microscopic origin for double-well poten- 
tials randomly distributed in Ge:0 systems; this is appar- 
ently in contrast with the phenomenological TLS model 
applied to amorphous solids, (ii) While the TLS model in 
amorphous solids requires an artificial distribution func- 
tion for two-level tunneling states, our model naturally 



yields the power-law distribution function with no artifi- 
cial condition. Consequently, the clarity of the coupled- 
rotor model enables us to argue the low-temperature 
properties of Ge:0 quantitatively, giving the theoretical 
prediction of the power-law specific heat in Ge:0. This is 
also the case for alkali haride crystals containing tunnel- 
ing dipoles. 



VII. CONCLUSION 

In conclusion, we have theoretically investigated the ef- 
fect of the dipolar interaction of oxygen defects on various 
physical properties in crystalline germanium. Numerical 
simulations have revealed that the dipolar interaction of 
nearby oxygen impurities engenders non-trivial low-lying 
excitations, causing anomalous behaviors for Ge:0 sys- 
tems with an oxygen concentration 10 17-18 cm~ 3 . The 
following were the main findings for the low-temperature 
anomalies of Ge:0: (i) the Schottky peak in specific heats 
C(T) at around 1 K, (ii) the power-law temperature de- 
pendence below O.IK, (iii) the p 2 -dependence of the mag- 
nitude of C{T) below 0.1 K, (iv) the anomalous hump 
in the dielectric susceptibility x(T) at around 1 K, and 
(v) the p 2 -dependence of the magnitude of the devia- 
tion A\(T) from the susceptibility of isolated rotors. All 
these behaviors originate from the presence of interacting 
quantum rotors in Ge:0 systems, and can be understood 
quantitatively by considering the contribution of the en- 
semble of interacting rotors. We have accounted for the 
occurrence of non-trivial low-lying excitations in inter- 
acting rotors in line with the two-level tunneling theory. 
The picture is based on quantum tunneling in a local 
double-well potential, thus making clear the relevance 
of the present results to other systems involving tunnel- 
ing constituents. We hope that our findings shed light 
on experimental research regarding the low-temperature 
properties of oxygen-doped semiconductors as well as on 
disordered systems with interacting dipole moments. 



Acknowledgments 

This work was supported in part by a Grant-in-Aid 
for Scientific Research from the Japan Ministry of Edu- 
cation, Science, Sports and Culture. One of the authors 
(H.S) thanks the financial supports from the Foundation 
Advanced Technology Institute. A part of numerical cal- 
culations were performed on the SR8000 of the Super- 
computer Center, Institute of Solid State Physics, Uni- 
versity of Tokyo. 



1 W. Hayes and D. R. Bosomworth, Phys. Rev. Lett. 23, and G. D. Watkins, Proc. R. Soc. A, 317, (1970) 133. 

(1969) 851; D. R. Bosomworth, W. Hayes, A. R. L. Spray 



13 



2 H. Yamada-Kaneta, C. Kaneta and T. Ogawa, 
Phys. Rev. B 42, 9650 (1990). 

3 F. Shimura, Oxygen in Silicon, Semiconductors and 
Semimetals, Vol. 42 (Academic Press, New York, 1994) 
and references therein. 

4 J. Kato, K. M. Itoh, H. Yamada-Kaneta and H-J. Pohl, 
Phys. Rev. B 68, 035205 (2003). 

5 W. Kaiser, J. Phys. Chem. Solids 23, 255 (1962). 

6 J. W. Corbet, R. S. McDonald, and G. D. Watkins, 
J. Phys. Chem. Solids 25, 873 (1964). 

7 P. Clauws, Mater. Sci. Eng., B 36, 213 (1996). 

8 M. Pesola, J. von Boehm, V. Sammalkorpi, T. Mattila, 
and R. M. Nieminen, Phys. Rev. B 60, R16267 (1999). 

9 H. Ch. Alt, Appl. Phys. Lett. 54, 1445 (1989); 55, 2736 
(1989); Phys. Rev. Lett. 65, 3421 (1990). 

10 J. Schneider, B. Dischler, H. Seelewind, P. M. Mooney, 
J. Lagowski, M. Matsui, D. R. Beard, and R. C. Newman, 
Appl. Phys. Lett. 54, 1442 (1989). 

11 M. Gienger, M. Glaser and K. Lassmann, Solid State Com- 
mun. 86, 285 (1993). 

12 A. J. Mayur, M. D. Sciacca, M. K. Udo, A. K. Ramdas, 
K. Itoh, J. Wolk, and E. E. Haller, Phys. Rev. B 49, 16293 
(1994). 

13 E. Artacho, F. Yndurain, B. Pajot, R. Rami'rez, C. P. Her- 
rero, L. I. Khirunenko, K. M. Itoh, and E. E. Haller, Phys. 
Rev. B 56, 3820 (1997). 

14 B. Pajot, P. Clauws, J. L. Lindstrom, and E. Artacho, 
Phys. Rev. B 62, 10165 (2000). 

15 J. Coutinho, R. Jones, P. R. Briddon and S. Oberg, Phys. 
Rev. B 62, 10824 (2000). 

16 H. Shima and T. Nakayama, Phys. Rev. B 69, 035202 
(2004). 

17 H. Shima and T. Nakayama, Phys. Rev. A 70, 013401 
(2004). 

18 See, Tunneling Systems in Amrphous and Crystalline 
Solids, edited by P. Esquinazi (Springer- Verlag, Berlin 
1998). 

19 An exact value of \p\ may be obtained through a quantum- 
chemical calculation; to perform this is beyond the present 
study. The influence of the change in the definition of \p\ to 
the results reached in this work will be presented at need. 

20 E. T. Whittaker and G. N. Watson, A course of modern 
analysis, (Cambridge University Press 1962). 

21 R. F. Gloden, Euratom Rpt. EUR 4349f (1970). 

22 See, A. Wiirger, Coherent Tunneling to Relaxation, 
(Springer- Verlag, 1997). 

23 M. P. Solf and M. W. Klein, Phys. Rev. B 46, 8147 (1992); 
47, 11097 (1993). 



24 P. W. Anderson, B. I. Halperin, and C. M. Varma, Phil. 
Mag. 25, 1 (1972). 

25 W. A. Phillips, J. Low. Temp. Phys. 7, 351 (1972). 

26 L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 
(Pergamon Press Oxford, 1977). 

27 W. Schnelle and E. Gmelin, J. Phys. Condens. Matter 13, 
6087 (2001). 

28 See, M. L. Mehta, Random Matrices (Academic Press, NY, 
1991), for derivating the explicit form of P(R) in one- 
dimensional systems. 

29 M. W. Klein, Phys. Rev. B 29, 5825 (1984); 31, 1114 
(1985). 

30 C. Enss, M. Gaukler, S. Hunklinger, M. Tornow, R. Weis, 
and A. Wiirger, Phys. Rev. B 53, 12094 (1996). 

31 A. Wiirger, R. Weis, M. Gaukler and C. Enss, Euro- 
phys. Lett. 33, 533 (1996). 

32 A. Wiirger, Z. Phys. B 94, 173 (1994); 98, 561 (1995). 

33 O. Terzidis and A. Wiirger, Z. Phys. B 94, 341 (1994). 

34 In our previous paper 35 , the contribution of paired rotors 
to x(T) was overestimated; this is due to that all rotors 
contained in Ge:0 were factorized into pairs regardless of 
a random distribution of rotors in a cubic space. 

35 H. Shima and T. Nakayama, J. Phys. Soc. Jpn 73, 2464 
(2004). 

36 S. D. Baranovskii, B. I. Shklovskii, and A. L. Efros, 
Sov. Phys. JETP 51, 199 (1980). 

37 A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975); 
in Electron-Electron Interaction in Disordered Systems, 
edited by A. L. Efros and M. Pollak, Elsevier Science Pub- 
lications B.V. (1985). 

38 S. K. Kirkpatrick and C. M. Varma, Solid State Commun. 
25, 821 (1978). 

39 R. Fisch, Phys. Rev. B 22, 3459 (1980). 

40 R. C. Zeller and R. O. Pohl, Phys. Rev. B 4, 2029 (1971). 

41 D. J. Salvino, S. Rogge, B. Tigner and D. D. Osheroff, 
Phys. Rev. Lett. 73, 268 (1994). 

42 A. L. Burin, J. Low. Temp. Phys. 100, 309 (1995); 
A. L. Burin and Y. Kagan, JETP 80, 761 (1995). 

43 D. Rosenberg, P. Nalbach and D. D. Osheroff, 
Phys. Rev. Lett. 90, 195501 (2003); S. Ludwig, P. Nal- 
bach, D. Rosenberg, and D. Osheroff, Phys. Rev. Lett. 
90, 105501 (2003); S. Ludwig and D. D. Osheroff, 
Phys. Rev. Lett. 91, 105501 (2003). 

44 J. Le Cochec and F. Ladieu, Eur. Phys. J. B 32, 13 (2003); 
F. Ladieu, J. Le Cochec, P. Pari, P. Trouslard, and P. Ail- 
loud Phys. Rev. Lett. 90, 205501 (2003). 



